function [centers, radius] = findIris(I)
close all;
    Image=I;
    G=rgb2gray(I);
    h=fspecial('gaussian',10,5);
    G1=imfilter(G,h);

    BW=edge(G1,'canny');

    figure(1);imshow(G1);
    [centers,radii,metric]=imfindcircles(BW,[20,60]);
    
    figure(2);
    imshow(I);hold on;
    viscircles(centers,radii,'EdgeColor','b');
    
    figure(3);
    imshow(BW);
end